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Abstract 



We estimate the variance and the skewness of the smoothed density field of simulated clusters 
in several dark matter (DM) models. The cluster simulations are based on the Zel'dovich 
approximation, the low computational cost of which allows us to run as many as 50 random 
realizations of each model. The simulated cluster distributions are smoothed with a Gaussian 
window at two different smoothing scales, R sm = 20 and 30/i~ 1 Mpc. We compare our results 
with a similar analysis of a redshift sample of Abell/ACO clusters (Plionis & Valdarnini 
1994). Within the list of considered models, we find that only the mixed Cold+Hot DM 
model (with fihot = 0.3) provides a good fit to the data. The standard CDM model and 
low-density (f2 = 0.2) CDM models, both with and without a cosmological constant term 
(Qa = 0.8), are ruled out. The tilted CDM model with primordial spectral index n = 0.7 
and a low Hubble constant [h = 0.3) CDM model are only marginally consistent with the 
data. 



Subject headings: Cosmology: large-scale structure of the universe - galaxies: clusters: 
general 



1 



1 Introduction 



On scales traced by galaxy clusters 20/i~ 1 Mpc), gravitational clustering has not yet 
entered its non-linear phase of evolution, so that the cluster distribution is mainly determined 
by the shape of the primordial fluctuation spectrum 1 . For this reason, several attempts 
have been devoted in compiling extended cluster redshift surveys both in the optical (e.g., 
Postman, Huchra k Geller 1992; Dalton et al. 1994, and references therein) and in the X-ray 
(e.g., Nichol, Briel k Henry 1994) bands. In order to compare the observational data sets 
with different cosmological models several authors have resorted to large N-body simulations 
which were designed to sample the length-scales relevant to the cluster distribution (e.g., 
White et al. 1987; Bahcall k Cen 1992; Croft k Efstathiou 1994). Since, however, linear, 
or at most mildly non-linear, gravitational clustering governs the cluster distribution, an 
N-body based approach is probably more sophisticated than what is really needed. On the 
other hand, existing analytical approaches, based either on linear theory (e.g., Bardeen et al. 
1986; Holtzman k Primack 1993) or on the Zel'dovich approximation (cf., Doroshkevich k 
Shandarin 1978; Mann, Heavens k Peacock 1993) rely on simplifying assumptions which do 
not allow an accurate treatment of the high-order correlation statistics or a realistic account 
to be taken of the observational biases, present in real data sets. 

In this Letter we analyze cluster simulations based on the Zel'dovich approximation 
(Zel'dovich 1970). It has been shown that this kind of simulation gives a reliable description 
of the cluster distribution at a fraction of the computational cost of N-body simulations 
(Borgani, Coles k Moscardini 1994; Borgani et al. 1994b; see also Blumenthal, Dekel k 
Primack 1988). Taking advantage of this, we analyze six different initial power-spectra and, 
at the same time, run a large number (50) of realizations of each model to get a reliable 
estimate of the cosmic variance on the correlation statistics. 

We analyze the density field traced by cluster distributions, after suitably smoothing it 
at different scales. This procedure has been recently applied to the cluster distribution by 
different authors (Plionis k Valdarnini 1994; Kolatt, Dekel k Primack 1994), since it limits 
the effect of shot-noise in clustering measures while a different approach, based on counts in 
cells and on measures of volume-averaged correlation functions, has been applied by Cappi 
k Maurogordato (1994). We evaluate the variance and the skewness of the resulting cluster 
density field. The comparison between our results and those obtained by Plionis k Valdarnini 
(1994; PV94 hereafter), using a similar analysis on an Abell/ACO redshift sample, allows us 
to put stringent constraints on the dark matter (DM) models we have considered. This is 
a further confirmation (see Borgani, Coles k Moscardini 1994) that our method to generate 
realistic cluster distributions represents a flexible instrument which can be used effectively to 
test the parameter space of DM models against the observed large-scale cluster distribution. 

1 h is the Hubble constant H in units of 100 km s _1 Mpc -1 . 
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The layout of this Letter is as follows: in Section 2 we describe briefly our simulations and 
the considered models; in Section 3 we present the analysis method and the results obtained; 
we state our main conclusions in Section 4. 



2 The Simulations 

Here we simply sketch our simulation procedure, based on the Zel'dovich approximation [ZA] 
(cf. Shandarin & Zel'dovich 1989); more details are given in Borgani et al. (1994b). 

Let q and x(q, t) be the initial (Lagrangian) and the final (Eulerian) comoving position 
of a fluid element. The ZA assumes that, at the time t, the Eulerian position is given by 



x(q,*)= [q+&(*) Vq^(q) 
perturbations and V'(q) is t 



Here b(t) is the growing mode of the evolution of linear density 
the gravitational potential, which is related to the initial (linear) 
density fluctuation field, £(q), through the Poisson equation. Accordingly, the fluid particles 
move under the ZA along straight lines. Although gravity determines the initial "kick" 
given to the fluid particles, they do not afterwards feel any tidal effects. Particles fall into 
gravitational potential wells to form structures, but these structures evaporate subsequently 
since no account is taken of the self-gravity of these non-linear structures. In reality particles 
will attract each other when their orbits cross or they pass close to each other but, in the ZA, 
matter always moves along its initial path, oblivious to the presence of the other particles. 
In this sense, the ZA gives a good description of gravitational dynamics as far as particle 
trajectories do not intersect with each other, i.e. before "shell-crossing": after this moment, 
the validity of the ZA breaks down. 

In order to improve the performance of ZA, one can filter out the small-scale wavelength 
modes in the linear power-spectrum, P(k), which are responsible for shell-crossing (Coles, 
Melott & Shandarin 1993; see also Kofman et al. 1992). Accordingly, we convolve the 
linear power-spectrum with a Gaussian filter, P(k) — > P(k) e~^ kR ^ 2 , which has been shown 
(Melott, Pellman & Shandarin 1994) to be the optimal choice. Kofman et al. (1994) provided 
an analytical expression for the average number of streams, N s , at each Eulerian point as 
a function of the r.m.s. fluctuation amplitude, cr, of the Gaussian density field. The limit 
corresponding to the single-stream regime, N s — > 1, is attained for a — > 0, with N s increasing 
rapidly at larger values of a: see eq.(27) in Kofman et al. (1994). We choose the filtering 
radius Rf for each model so that N s = 1.1. We found this choice to be a reasonable 
compromise between avoiding shell-crossing and suppressing the development of genuine 
clustering. In the second column of Table 1 we report the filtering radii used for the DM 
models considered in this paper. 

We have generated random-phase realizations of the linear density field on 128 3 grid 
points for a box of side L = 320/i~ 1 Mpc. We moved particles from initial grid positions 
according to the ZA and re-assigned the density and the velocity field on the grid through 
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a TSC interpolation scheme (see e.g. Hockney & Eastwood 1981) for the mass and the 
momentum carried by each particle. Then we defined clusters as local density maxima on 
the grid according to the following prescription. If d c i is the average cluster separation, 
then we select N c i = (L/d c i) 3 clusters at the N c i highest density peaks. In the following, 
we assume d c i = 40/i~ 1 Mpc, which is the corresponding value of the combined Abell/ACO 
cluster sample to which we will compare our simulation results (see PV94). 

We ran simulations for six different models of the initial fluctuation spectrum. For each 
model we generated 50 random realizations, in order to estimate the cosmic variance reliably. 
The models we have considered are the following: the standard CDM model (SCDM) with 
cr 8 = 1 for the r.m.s. fluctuation amplitude within a top-hat sphere of 8/i _1 Mpc; a tilted 
CDM model (TCDM), with primordial spectral index n = 0.7, h = 0.5 and cr 8 = 0.5 (e.g., 
Cen et al. 1992; Tormen et al. 1993; Adams et al. 1993; Moscardini et al. 1994); a 
low Hubble constant CDM model (LOWH), with h = 0.3 and a 8 = (1.6) -1 (Bartlett et 
al. 1994); a Cold+Hot DM model (CHDM), with fi hot = 0.3 for the fractional density of 
the hot component, h = 0.5 and cr 8 = (1.5) -1 (e.g., Klypin et al. 1993); an open CDM 
model (OCDM), with Q = 0.2, h = 1 and c 8 = 1; a spatially flat, low-density CDM model 
(ACDM), with Q = 0.2, h = 1, a 8 = 1.3 and cosmological constant = 0.8 (Bahcall & Cen 
1992; Dalton et al. 1994). Transfer functions have been taken from Holtzman (1989), except 
that of LOWH, which was taken from Bond & Efstathiou (1984), with suitably chosen shape 
parameter T = Q h = 0.3. All the models, except the open CDM one, are normalized to 
be consistent with the quadrupole of the CMB temperature anisotropy measured by COBE 
(Bennett et al. 1994). In a forthcoming paper (Borgani et al. 1994b), we will present more 
details about our simulations. 

3 The analysis 

In the following, we implement the same analysis as in PV94, in order to allow a reliable 
comparison between real cluster data and our simulations (we refer to that paper for further 
details on the the real cluster sample and the analysis method; see also Borgani et al. 1994b). 
We obtain a continuous cluster density field by smoothing the cluster distribution on a grid, 
with grid-cell width of 20 h~ x Mpc, using a Gaussian kernel. In order to study the cluster 
density field at different smoothing scales, i? sm , we use two radii for the Gaussian kernel, 
namely R sm = 20 and 30 h~ x Mpc. If p(x g ,R sm ) is the smoothed density field at the grid 
point position x g , then the relative fluctuations are given by £(x g , R sm ) = p(x g ,R sm )/p — 1, 
where the average density p does not depend on R srn . The variance a 2 (R sm ) and the skewness 
~f(R S m) are defined as the second- and third-order moments of the 6 field, respectively. We 
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therefore have 

a\R g ) = ^-J2l6(x g ,R sm )f ; 7W = ^E[^)f. (!) 

Iy gr g =i ly gr g= i 

where N gr = (16) 3 is the total number of grid points. We have chosen not to take into 
account Poisson shot-noise corrections since the cluster distribution can be hardly considered 
a Poisson sampling of the underlying (galaxy) distribution (cf. Coles & Frenk 1991; Borgani 
et al. 1994a). Moreover, Gaztahaga & Yokoyama (1993) have shown that the smoothing 
process itself suppresses considerably the shot-noise effects. For these reasons PV94 did not 
use any shot-noise corrections so, to make a consistent comparison of our models with their 
results, we did not include such corrections in our analysis either. Furthermore, since all 
the model cluster distributions have the same (p) and we treat them similarly, the possible 
effects of shot-noise are relatively and qualitatively canceled out in the model and data 
intercomparison. 

In Figure 1 we plot the results of our analysis in the a 2 -^ plane, following Coles & Frenk 
(1991; see also Coles et al. 1993). In each panel, we show the variance-skewness relation 
estimated in redshift space and at the two smoothing scales for the corresponding model. 
The analysis has been done in redshift space in order to compare our results properly with 
real data. Redshift distortions affect the a 2 and 7 values by ~ 10% at R sm = 20/i~ 1 Mpc 
and by ~ 6% at R sm = 30/i~ 1 Mpc (see also Borgani et al. 1994b). Each point represents 
one of the 50 realizations, so that the scatter represents the effect of cosmic variance. The 
crosses represent the results of the PV94 analysis of the combined Abell/ACO cluster redshift 
sample. In Table 1 we report, for each model, the values of a 2 and 7 estimated as the 
average between the 50 realizations and their "cosmic r.m.s." uncertainties; we also give the 
corresponding PV94 results. 

From Figure 1 we see that R sm = 20/i~ 1 Mpc is the most effective of the two scales we 
have considered for discriminating between the models. At this scale only the CHDM model 
fits best the data. The LOWH and TCDM are marginally consistent, with only about one 
tenth of the realizations giving a 2 and 7 values as high as those of the Abell/ACO sample. 
The SCDM, ACDM and OCDM are ruled out, giving rise to either a too weak (SCDM) 
or a too strong clustering (OCDM and ACDM). The results at R sm = 30/i~ 1 Mpc are less 
discriminative. At this scale, the only models which can be ruled out to a high confidence 
level are the ACDM and OCDM for which, again, no realization reproduce the PV94 result. 
Although the other models are all consistent with observations, CHDM and TCDM seem to 
fare better (see also Table 1). We can conclude, therefore, that R sm = 20/i~ 1 Mpc is a sort of 
optimal scale where to test dark matter models against the distribution of galaxy clusters. 
At smaller scales shot-noise dominates in the estimation of the cell-count moments, while 
at larger scales the clustering becomes significantly weaker, and consequently more difficult 
to measure from a finite sample, be it real or simulated. 
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4 Discussion & Conclusions 



We have used extensive simulations of rich galaxy clusters, based on the Zel'dovich approx- 
imation. We considered six different models for the linear power-spectrum and generated 
a large number (50) of random-phase realizations of each model, in order to determine ac- 
curately the effect of cosmic variance. We then generated a smoothed density field on a 
grid by convolving the cluster distributions with a Gaussian window at two different scales, 
R sm = 20 and 30/i~ 1 Mpc. In order to constrain the power-spectrum models, we computed 
the variance and the skewness of the density field and compared it with the results of a 
similar analysis of an Abell/ACO combined cluster redshift sample, realised by Plionis & 
Valdarnini (1994). 

We find that R sm = 20 h~ 1 M.pc is the most reliable scale at which to discriminate between 
models. At this scale, only the CHDM model fits the data very well. SCDM, ACDM and 
OCDM are ruled out at a high significance level. The LOWH and TCDM models are 
marginally consistent with the data. At R sm = 30/i~ 1 Mpc the only models to be definitely 
ruled out are ACDM and OCDM. Since the OCDM and ACDM have a very similar P(k) 
shape, the higher normalization of the latter is the reason for the even stronger clustering it 
produces with respect to the former. 

Clearly, one can imagine making modifications to the choice of the parameter values, 
involved in the definition of the DM models we have considered, in such a way as to make 
them consistent with the data. The physical justification for such an ad hoc tuning of 
parameters is, however, quite feeble at the present time. For example, in the LOWH model, 
it is quite hard to justify at the present time a Hubble parameter as small as h = 0.3 
(Fukugita, Hogan & Peebles 1993). As far as TCDM is concerned, lowering the primordial 
spectral index still further, below n ~ 0.7, could increase the excess power on large scales. 
However, a value even as low as n = 0.7 is strongly constrained by the level of CMB 
anisotropy detected by COBE (Bennett et al. 1994) and by the properties of the large-scale 
galaxy peculiar velocity field (e.g., Moscardini et al. 1994). Increasing the density parameter 
in the OCDM and ACDM models should decrease the cluster correlation to a satisfactory 
level: we have, in fact, verified that simulations based on Q = 0.4 fare much better than 
those we have considered here. In a forthcoming paper, we plan to exploit appropriate 
properties of the cluster velocity field in order to discriminate even more between high- 
density and low-density universes. 

As a final comment we should mention that, although the CHDM model with fihot = 0.3 
is consistent with cluster correlation data, it has nevertheless been shown to be severely con- 
strained by the detection of high-redshift objects (e.g. Klypin et al. 1994; Ma & Bertschinger 
1994). Although lowering the hot fraction to fihot = 0.2-0.25 may alleviate problems con- 
nected with the timing of galaxy formation in this scenario, it remains to be seen whether 
this can be achieved without destroying the apparent consistency of the model with the 
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observed clustering of rich clusters that we have found in this work. 

In the light of these further possible investigations, it is clear that the ease and low 
computational cost of our method for generating cluster simulations, provides an extremely 
useful tool for investigating the large parameter space of dark matter models and comparing 
them with real data. 
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Table 1: Column 2: Filtering radius (in units of h~ 1 M.pc for the linear power-spectrum. 
Column 4 and 5: Variance and skewness in redshift space, at the two smoothing radii 
indicated in column 3 (in units of /j, _1 Mpc), for the six simulated models and for the real 
data (from PV94). 



Model 



Rf R s 
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SCDM 4.4 20 0.373±0.029 0.262±0.060 

30 0.110±0.014 0.023±0.011 

TCDM 1.6 20 0.385±0.044 0.297±0.097 

30 0.123±0.024 0.032±0.018 

LOWH 2.4 20 0.385±0.042 0.282±0.080 

30 0.119±0.020 0.027±0.014 

CHDM 2.2 20 0.504±0.050 0.499±0.146 

30 0.155±0.024 0.050±0.027 

OCDM 4.5 20 0.638±0.058 0.771±0.189 

30 0.215±0.031 0.089±0.039 

ACDM 6.3 20 0.899±0.095 1.689±0.689 

30 0.298±0.057 0.211±0.168 



Abell/ACO 
(PV94) 



20 0.496±0.064 0.445±0.187 
30 0.135±0.028 0.033±0.028 
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Figure Caption 

Figure 1. Scatter-plots for the variance-skewness relation. Each point refers to a single 
realization of the corresponding model. Triangles and squares are for R sm = 20/i~ 1 Mpc and 
Rsm = 30/i _1 Mpc, respectively. Heavy crosses are the results of a similar analysis realized 
by PV94 on a combined Abell/ACO cluster redshift sample. 
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